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Atomic bosons and fermions in an optical lattice can realize a variety of interesting condensed 
matter states that support equilibrium current patterns in the presence of synthetic magnetic fields 
or non-abelian gauge fields. As a route to probing such mass currents, we propose a nonequilibrium 
quantum quench of the Hamiltonian that dynamically converts the current patterns into experi- 
mentally measurable real-space density patterns. We illustrate how a specific such "unidirectional" 
quench of the optical lattice can be used to uncover bulk checkerboard and stripe current orders in 
lattice Bose superfiuids and Fermi g well as chiral edge currents in a quantum Hall state. 



Quantum dynamics of particles moving in magnetic 
fields or non-abelian gauge fields is of great interest in 
various condensed matter systems such as the cuprate 
superconductors [T] , quantum Hall liquids [2] , and topo- 
logical insulators [2j ■ Such gauge fields can result in equi- 
librium charge or spin currents of electrons. For instance, 
a type-II superconductor in a magnetic field forms an 
Abrikosov vortex lattice which supports a periodic bulk 
current pattern [1]. A uniform magnetic field for lattice 
electrons can lead to topologically nontrivial states with 
a quantized Hall conductance and chiral edge currents 
[5]. In a solid, such electronic currents produce their own 
characteristic magnetic fields, and can thus be probed by 
using magnetic microscopy [6] or neutron scattering 7\. 

Recently, experiments have begun exploring effects of 
"artificial" orbital magnetic fields [8HTT| and non-abelian 
gauge fields [T^ on neutral ultracold atomic gases. These 
experiments can potentially realize a wide variety of 
fermionic and bosonic states with equilibrium mass cur- 
rents. This brings us to a crucial question of great in- 
terest: How can we observe the equilibrium mass current 
patterns for neutral atomic gases? 

In this paper, we propose a dynamical probe of cur- 
rents, which relies on experimental progress in measur- 
ing lattice scale modulations of the atom density. Such 
density mapping tools include noise correlations [13l [M] , 
Bragg scattering |15j . and in situ microscopy at the lat- 
tice scale [TB]. Our key idea is to make a specific quan- 
tum quench of the Hamiltonian that violates the steady 
state divergence-free condition on equilibrium currents, 
thus causing an imbalance between currents entering and 
leaving different sites. This leads to characteristic density 
build-up or depletion, as dictated by the continuity equa- 
tion. Imaging the subsequent density variation across the 
lattice then yields real space information about the initial 
currents. Our proposal thus links up with another very 
important thread in ultracold atom research: quench- 
induced quantum dynamics in many-body systems |17| . 

Consider a thought experiment in which we start in 
an equilibrium state with currents and instantaneously 
turn off the particle hopping amplitude on every bond 
except one. Then, the instantaneous current J7 on that 
one bond is unaffected, but currents on all other bonds 



vanish. A short time St after this local quench, there will 
be an accumulated density imbalance J^St between the 
two sites connected by the unperturbed bond; measuring 
this yields information on the magnitude and direction 
of the bond current. A similar such "quasi- local" current 
probe has been used in a recent experimental study of 
nonequilibrium dynamics in a ID Bose gas |18j . 

In this paper, we study equilibrium current patterns 
in 2D, focusing on a specific "unidirectional" quench, 
where we suddenly decrease the particle hopping am- 
plitude in one direction across the entire lattice. This 
is achieved simply by changing the optical lattice laser 
intensity along one direction. We show that tracking 
the subsequent space-time variation of the density for 
this quench protocol can uncover the underlying current 
patterns in many interesting cases which have been cho- 
sen to show that the method works for both bosons and 
fermions, and can probe bulk as well as edge currents. 
Specifically, we illustrate this scheme for (i) a square lat- 
tice Bose superfluid in a staggered (checkerboard) mag- 
netic flux pattern, (ii) a square lattice Bose superfluid 
with a striped magnetic flux pattern, (iii) noninteract- 
ing fermions on a square lattice in a staggered magnetic 
flux, and (iv) the integer quantum Hall state of lattice 
fermions in a uniform magnetic fiux that supports chiral 
edge currents. Such quenches can also probe the current 
pattern in the recently studied chiral Bose Mott insulator 
119' and spin currents of atomic matter. 

Bosons in a checkerboard flux pattern. — Consider in- 
teracting bosons on a 2D square lattice, described by 
the Bose Hubbard model Hbh — ~J2rr''^rr'^lK' + 

X^rK ^r^r + ^ X^r ^r^r^r^r ■ Here Jr^r' dcuOtCS the bo- 

son hopping amplitude between sites (r,r'), U is the 
on-site repulsion, and Vr is a harmonic trap potential 
given by K = Vo{x^ + V^)- We take J^y = J^, ^ ^ 
only for nearest neighbors, and choose Jr,r+s — J and 
Jr,r+y = Jy exp(z(— 1)^+''(/)/2) . This yields staggered 
magnetic fluxes, ±(j), that pierce the elementary square 
plaquettes in a checkerboard pattern; a route to realizing 
such a flux pattern has been proposed previously pi] . 

For weak interaction, U <J,Jy, we first solve for the 
equilibrium ground state |22j for Jy = J by numerically 
minimizing the Gross-Pitaevskii (GP) energy fimctional 
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m EgP = -Er,r' -/r^r'^J^r' +Er(|l*r|^ + l^l^rP), 

where is the condensate wavefunction at lattice site 
r = {x, y). As shown in FigQa) for a system with hnear 
length L = 22, Vq = 0.05J, and an average filling factor 
of n — 1, this leads to a superfluid ground state with 
staggered loop currents. The smooth density profile of 
the ground state reflects the trap potential, but it does 
not reveal the currents induced by the gauge field. 
Starting with this ground state, we suddenly decrease 



Jy from its initial value J* 



J to a final value < J 



at time t = 0. The subsequent condensate dynamics is 
described by the time-dependent GP equation pU] 



ih 



dt 



J2j!,r''^r'{t) + {Umt)WK)Mt). (1) 



Henceforth, we set H = 1. To study the dynamics, we 
numerically integrate this time-dependent GP equation 
[22] , which allows us to extract the condensate wavefunc- 
tion ^r{t), and the boson density |^'r(t)P, at later times. 
As shown in Figjl] the condensate at t > develops a 
striking checkerboard density pattern that reflects the 
underlying current order. Information about the direc- 
tion of circulation on a plaquette is easily discerned from 
the density pattern established after a short time period; 
since the quench is along Jy , the initial build up of density 
is on sites that have currents flowing into them along the 
strong J-bonds oriented along the x-axis. After a short 
time has passed, the density buildup reaches a maximum 
and the flow is reversed, resulting in "plasma oscillations" 
between the two checkerboard patterns. The frequency of 
these oscillations scales as ~ y/UJ\W^-, it varies slowly 
with position due to the density inhomogeneity in the 
trap. Such checkerboard oscillations persist for a long 
time, as inferred from the density difference between the 
two sublattices AnAsit) = nA{t) —nsit) shown in Figjl] 

We find that the sublattice density pattern oscillations 
are also robust against moderate random phase fluctua- 
tions imprinted on the initial state [2^. This suggests 
that thermal phase fluctuations will suppress but not 
completely destroy these oscillations. Remarkably, the 
sublattice density oscillation persists out to long times, 
^ 1, as seen in Figjljb). At these times, shown 
in Fig. [l]^e) and (f), we find additional long- wavelength 
modulations superimposed on the checkerboard density 
pattern, (i) Spherical density waves emanate periodically 
outward from the centre, which we attribute to the spa- 
tial variation of the "plasma frequency" resulting from 
the radial variation of the density l^'rP in the trap, (ii) 
The cloud shape shows oscillatory distortions into an el- 
lipse due to the anisotropy of the final tunneling Jy < J. 

Bosons in a stripe flux pattern. — We next con- 
sider the above Bose Hubbard model in the presence 
of a striped magnetic flux pattern as realized in a re- 
cent experiment [5] . We choose Jr.r+x — J and Jr.r+y = 
Jyexp{i{—iy(j)x), so that we enclose fluxes ±0 through 
each plaquette that lies along a stripe in the vertical di- 
rection, and solve for the equilibrium ground state by 
minimizing the GP energy functional for U'^J,Jy, with 
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FIG. 1: (Color online) Density pattern of 2D Bose super- 
fluid in a staggered flux gauge field following a quench, (a) 
Initial density profile and (inset) current pattern of conden- 
sate ground state at half-filling, (b) Time dependence of 
the sublattice density difference An{t) for a staggered fiux 
of |(^| = 7r/2 flowing quench oi Jy = J to = 0.25J (dashed 
blue) and = 0.50 J (solid red), (c - f) Change in local 
density, 5n (relative to original density), at different times 
following a quench Jy = J ^ J^ = 0.50J with U = J. Cir- 
cled region in (c), (d) indicates central plaquette vortex. 



L = 22, Vq = 0.05J, and an average filling n — 1. We 
find a superfiuid with vertically striped loop currents de- 
picted in Figj2j which resembles a stripe pattern of 'long 
vortices' that are highly elongated along the y-direction. 
Again, the smooth equilibrium density pattern reveals no 
information about the underlying current order. 

Upon quenching Jy, the superfluid generates a den- 
sity pattern that strikingly reflects the underlying equi- 
librium striped currents. Each vertically elongated loop 
forms four quadrants of alternating high and low den- 
sity, giving rise to an oscillatory quadrupole moment 
SL°y ^ a;?/p(x, y, t) about the center of each 'long vortex', 
where the sum is restricted to a single vertically elongated 
loop composed of square plaquettes. The density pattern 
oscillates along the top (or bottom) half of the observed 
density profile as we go along the x-direction. At inter- 
mediate times, we again find that the spatial dependence 
of the oscillation frequency leads to additional long wave- 
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FIG. 2: (Color online) Dynamical density pattern of interact- 
ing bosons on a 2D square lattice in a stripe-like magnetic 
flux pattern following a quench, (a): Initial density pro- 
file and (inset) current pattern of condensate ground state 
at half-filling, (b - e): Change in local density, 5n (relative 
to the original density), at different times following a quench 
4 = J -> J/ = 0.50 J with U = J. Circled region indicates the 
central elongated vortex with a nonzero quadrupole moment. 



length stripe-like waves emanating from the trap centre. 
Quenching J (i.e., along the a;-direction) rather than Jy 
leads to a similar early time density patterns. However 
the oscillatory dynamics occurs on longer time scales due 
to mass transport occurring over larger distances ~ L. 

Spinless fermions in a staggered flux. — Motivated 
by exploring such quench-induced density dynamics for 
fermions, we next turn to noninteracting fermions in a 
staggered flux background [531. We study the Hamil- 

tonian Hsf = - I]r,r' •^r,r'/r /r'^ whcrC Jr.r+x = J 

and Jr,r+t/ — exp(i(— 1)'^+^(/)/2), leading to staggered 
checkerboard fluxes ±0 [5T]. In momentum space 
Hsf = EL^^k^Ucos^kT^ + sin0kT^)^'i^, where ^-j, = 

(/k' /k+o)' '''^'^ ^""^ Pauli matrices, and the prime on the 
sum implies that only momenta in the reduced Brillouin 
zone (BZ) are included. Here, we have defined 51k — 
Ek/i^k, and sin6'k = 7k/^^k, with ek = 



7k' cos( 



— 2(Jcos kx -\- Jy cos I cos ky) and 7k = — 2Jy sin | cos ky. 
This leads to mode energies irik- 

Imagine fermions initially filled into negative energy 
states — fik up to a Fermi energy Ep, and then quench- 
ing Jy from Jy—i'J^ at time t — 0, which instantaneously 

changes (i^ki^'k) (^^k,^k)- This translationally invari- 
ant quench ensures that different momentum pairs (k, k-|- 
Q) stay decoupled from each other. The ensuing "spin 
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FIG. 3: (Color online) Time dependence of the sublattice den- 
sity difference AnAs(t), for noninteracting spinless fermions 
on a 2D square lattice at a filling of n = 0.4, and various 
staggered fiuxes 0, following a quench Jy — J ^ 



Jl = 0.5J. 



precession" type dynamics in momentum space leads to 
an oscillating density difference between the two sublat- 
tices AnAB{t) = ^ E{k} J sin(6'k - ^k) sin(2f2ki) 
where the momentum sum runs over only initially occu- 
pied states in the reduced BZ. A numerical evaluation 
of the sum allows us to plot the sublattice density oscil- 



lations, shown in Fig. 3 for J, 



J , Jy 



0.5 J, fermion 



density n = 0.4 per site, and various staggered flux values. 

These oscillations exhibit multiple frequencies due to 
the large number of occupied fermion modes; the very 
generation of such a sublattice imbalance is evidence of 
an initial staggered current pattern. Over the entire 
range of displayed fluxes, and a wide range of densi- 
ties n ~ 0.3-0.5 near half-filling, we find that the dom- 
inant oscillation frequency arises from initially occupied 
states near k — (0, tt) due to a van Hove singularity 
in the density of states [52]. This leads to an esti- 
mated, nearly density-independent, dominant oscillation 

frequency £7* 



2^^ J2 + (4)2 - 2JJI cos f , in good 
agreement with the numerical data in Fig. [3] The weak 
density dependence of An^B(t) over a range of fillings 
indicates that trap induced inhomogeneities will not sig- 
nificantly affect these oscillations. 

Lattice quantum Hall state. — Finally, we turn to 
topologically nontrivial states that support edge currents. 
Recent work has focused on extracting the nontrivial 
band topology from time-of-fiight measurements [53 [5S] 
or spectroscopy of the edge modes [27]. Here we explore 
density dynamics induced by the unidirectional quench 
for lattice fermions in a uniform magnetic field. Propos- 
als to obtain such uniform fluxes exist in the literature 
[l8] [28] . For concreteness, consider fermions on a 2D 
square lattice with a uniform magnetic flux (/)=27r/3 per 
plaquette. The resulting particle-hole symmetric Hofs- 
tadter spectrum |29j has three non-overlapping bands, 
with Chern numbers -1-1, —2, -1-1, so that 'band insula- 
tors' with some bands being completely flUed support 
a nonzero quantized Hall conductance, and chiral edge 
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FIG. 4: (Color online) Density dynamics for spinless fermions 
in the lowest Hofstadter band with flux (j> = 2ii/'i per plaque- 
tte on a square lattice following a quench from Jy — J to 
Jy = 0.5 J. (a) The scaled quadrupole moment Qxy{t)/L^ = 
ii/L^)'}2x ^?/<^^(2^i J/j for various system sizes versus the 
scaled time tJ/L. (b - d) Stripe-like density modulations (for 
L = 24) moving from the y-edges into the initially incom- 
pressible bulk at different indicated times. 



currents, yielding lattice quantum Hall (QH) states fS|. 

We begin by numerically diagonalizing the Hamilto- 
nian Hqu = - Sr.r'-^r.r'/^/r' "^1*^ Jr,r+x = J and 



./r,r r,r' r r' 

Jr.r+v — JyS^'^^, for (f> = 27r/3, with open boundary con- 
ditions on a L X L system, and fill up the lowest band 
(and some edge modes) to get a fermion filling n = 1/3. 
We find that the ground state bulk density is uniform 
(see Fig. |4] for i = 0) and supports edge currents con- 
fined to an "edge layer" of thickness '-^ 2 — 3 lattice sites. 
We next track the density dynamics following a quench 
from Jy = J to < J, which is easy to study once we 
compute the initial and final spectrum and eigenstates. 

Viewing the chiral edge currents as analogous to that 
arising from a 'vortex', we expect the quench to lead 
to quadrupolar density oscillations and current reversals, 
similar to what we found for the 'long vortex' in the stripe 
fiux superfluid. Inspired by recent work on the superfluid 
Hall effect of atomic bosons [30j, we study the behavior 



of the quadrupole moment Qxy{t) = y ^y^''^i^^y^^)' 
where the density deviation Sn is with respect to fi= 1/3. 
We find that Qxyii) indeed displays oscillatory sign re- 
versals and, as seen in Fig|4]ja), the data for various L 
collapse when plotted as Qxy{t)/L^ versus t/L. The t/L 
scaling shows that the oscillations occur due to trans- 
port across the system length L. A simple scaling argu- 
ment for an edge current induced oscillation shows that 
Qxy ^ L^, as we also see numerically. The numerical ob- 
servations are thus consistent with the quadrupolar os- 
cillations being driven by edge currents. For (f) — — 27r/3, 
Qxy{t) has the opposite sign. Taken together, these ob- 
servations provide strong evidence that the initial state is 
a nontrivial insulator which is incompressible in the bulk 
and supports chiral edge currents. 

We observe two additional effects, (i) There are rapid 
small amplitude oscillations of Qxy (t) , over a time scale 
^ superimposed on the long time sign reversals, 

(ii) The quench leads to the dynamic stripe-like density 
modulations (see Fig. |4]) which nucleate near the y-edges 
and propagate inward into the bulk from both directions. 
We attribute both observations to quench-induced 'band 
mixing'. Such mixing occurs across the lattice, but the 
asymmetry at the edge leads to the stripes originating 
from the edge. Such mixing with bands of opposite Chern 
number is also crucial for edge current reversals |22| . 

Discussion. — We have shown that quantum quenches 
can yield a probe of underlying current patterns of atoms 
in an optical lattice by converting them into measur- 
able real-space density orders. A possible concern re- 
garding such sudden quenches is that it may invalidate 
our tight-binding model since particles can get excited 
to very high bands of the periodic optical potential. To 
address this concern, the quench needs to be "adiabatic" 
on time-scales comparable to the inverse interband gap 
while being "sudden" on time scales governing intraband 
dynamics. This seems to be achievable; for instance, time 
dependent lattice modulation experiments on repulsive 
atoms in deep lattices can access regimes where they are 
consistent with theory on the one-band Hubbard model 
[5T] . Looking ahead, it would be extremely useful to 
explore quench protocols for distinguishing various topo- 
logical states that may be realized in future experiments. 
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